J ' 1 

TECHNICAL 

LIBRARY 


TECHNICAL  REPORT  ARLCB-TR-78014 


*p*Ao<>0  <>/3 


TEMPERATURES  AND  STRESSES  DUE  TO  QUENCHING 
OF  HOLLOW  CYLINDERS 


John  D.  Vasilakis 


August  1978 


US  ARMY  ARMAMENT  RESEARCH  AND  DEVELOPMENT  COMMAND 

LARGE  CALIBER  WEAPON  SYSTEMS  LABORATORY 
BENET  WEAPONS  LABORATORY 
WATERVLIET,  N.  Y.  12189 


A' 


AMCMS  No.  6111019A0011 
DA  Project  No.  1L161101A91A 
Pron  No.  1A825567GGM7 


t 


APPROVED  FOR  PUBLIC  RELEASE;  DISTRIBUTION  UNILIMTED 


DISCLAIMER 


The  findings  in  this  report  are  not  to  be  construed  as  an  official 
Department  of  the  Army  position  unless  so  designated  by  other  author- 
ized documents. 

The  use  of  trade  name(s)  and/or  manufacturer (s)  does  not  consti- 
tute an  official  indorsement  or  approval . 


DISPOSITION 

Destroy  this  report  when  it  is  no  longer  needed.  Do  not  return  it 


to  the  originator. 


SECURITY  CL  ASS1 FICATIOU  OF  THIS  PAGE  (Wh*n  Date  Entered) 


REPORT  DOCUMENTATION  PAGE 

READ  INSTRUCTIONS 

BEFORE  COMPLETING  FORM 

1.  REPORT  NUMBER 

ARLCB-TR- 78014 

2.  GOVT  ACCESSION  NO. 

3.  RECIPIENT’S  CATALOG  NUMBER 

«.  TITLE  (mnd  Subtftte) 

TEMPERATURES  AND  STRESSES  DUE  TO  QUENCHING  OF 
HOLLOW  CYLINDERS 

5.  TYPE  OF  REPORT  & PERIOD  COVERED 

6.  PERFORMING  ORG.  REPORT  NUMBER 

7.  AU  THORf a) 

John  D.  Vasilakis 

8.  CONTRACT  OR  GRANT  NUMBER^*) 

9.  PERFORMING  ORGANIZATION  NAME  AND  ADDRESS 

Benet  Weapons  Laboratory 

Watervliet  Arsenal,  Watervliet,  N.Y.  12189 
nRnAR-J.HR-TT. 

10.  PROGRAM  ELEMENT.  PROJECT,  TASK 
AREA  * WORK  UNIT  NUMBERS 

AMCMS  No.  6111019A0011 

DA  Project  No.  1L161101A91A 
Pron  No.  1A825567GGM7 

11.  CONTROLLING  OFFICE  NAME  AND  ADDRESS 

US  Army  Armament  Research  and  Development  Command 
Large  Caliber  Weapon  Systems  Laboratory 

Dover,  New  Jersey  0780!). 

12.  REPORT  DATE 

August  1978 

13.  NUMBER  OF  PAGES 

2A 

14.  MONITORING  AGENCY  NAME  & ADDRESS^!/  different  from  Controlling  Office) 

13.  SECURITY  CLASS,  (of  thlt  report) 

UNCLASSIFIED 

15*.  DECLASSI  FI  C ATI  ON/DOWN  GRADING 
SCHEDULE 


16.  DISTRIBUTION  STATEMENT  (of  thia  Report) 


Approved  for  public  release;  distribution  unlimited. 


17.  DISTRIBUTION  STATEMENT  (of  the  ebmtrect  entered  In  Stock  20,  It  different  from  Report) 


18.  SUPPLEMENTARY  NOTES 


19.  KEY  WORDS  (Continue  on  reverme  aide  if  necememry  end  Identify  by  block  number) 

Finite  Differences 
Transformation  Stresses 
Transient  Heat  Conduction 


20.  ABSTRACT  (Continue  on  reverem  aide  if  neceeaary  and  identify  by  block  number) 

After  forging,  gun  tube  blanks  are  heated  to  a high  temperature  and  quench- 
ed to  near  room  temperature  before  tempering  to  achieve  the  required  material 
properties.  The  purpose  of  the  quench  is  to  bypass  the  knee  of  the  pearlite 
phase.  This  program  was  undertaken  to  establish  cooling  curves  while  the  mater- 
ial is  being  quenched  and  to  compute  the  thermal  and  transformation  stresses 
involved. 

Continued  on  reverse  side 


DD  ,xr,»  1473  EDITION  OF  T NOV  65  IS  OBSOLETE 


SECURITY  CLASSIFICATION  OF  THIS  PAGE  (When  Date  Entered) 


SECURITY  CLASSIFICATION  OF  THIS  PAGE^on  Beta  Entered) 


Continued  from  Block  20. 

The  temperatures  are  computed  using  implicit  finite  difference  schemes. 

The  problem  treated  is  a nonlinear  one  in  radial  heat  flow.  The  problem  with 
cylindrical  geometry  is  assumed  to  be  axisymmetric  and  the  coefficients  in 
the  equation  such  as  thermal  conductivity  are  treated  as  functions  of  tempera- 
ture. The  boundary  conditions  are  written  in  a general  form  allowing  the  use 
of  temperature,  convection  or  heat  flux  boundary  conditions.  The  nonlinear 
problem  is  solved  by  using  two  finite  difference  schemes  in  tandem.  The  first 
computes  the  temperatures  at  the  n+h  time  step  assuming  constant  coefficients 
computed  from  a previous  temperature  distribution.  This  generates  a temperature 
distribution  throughout  the  thickness  which  is  used  to  compute  new  coefficients 
for^  the  second  finite  difference  scheme  which  calculates  the  temperature  dis- 
tribution at  the  n+1  time  step.  This  process  is  continued  until  a steady  state 
or  some  desired  level  is  reached. 

At  each  time  step,  the  program  computes  the  thermal  stresses  associated 
with  the  temperatures.  In  addition  to  this,  when  the  temperature  reaches  a 
certain  level,  called  martensite  start  (Ms) , the  material  begins  to  undergo  the 
martensite  transformation.  This  transformation  involves  an  increase  in  material 
volume  of  about  3%-4%.  A simple  view  of  these  transformation  stresses  is  taken 
and  the  stresses  due  to  this  volume  change  are  computed  as  the  temperature  cools 
to  below  the  martensite  start  temperature  throughout  the  wall  thickness. 

Results  are  presented  for  various  boundary  conditions  including  those 
expected  to  exist  in  the  quenching  facility. 


SECURITY  CLASSIFICATION  OF  THIS  P AGEfWTien  Data  Entered) 


ACKNOWLEDGEMENT 


The  author  wishes  to  acknowledge  the  assistance  of  Richard  Haggerty 
and  Royce  Soanes  in  the  computational  phase  of  this  work. 


i 


TABLE  OF  CONTENTS 


Page 

Acknowledgement  i 

I.  Introduction  1 

II.  Problem  Statement  2 

III.  Finite  Difference  Equations  3 

IV.  Thermal  and  Transformation  Stresses  7 

V.  Results  and  Discussion  9 

References  10 

LIST  OF  ILLUSTRATIONS 

Figure 

1.  Node  Placement  with  Points  Shifted  from  Boundaries  11 

2.  Typical  Temperature  Distributions  Near  Transformation 

Temperature  12 

3.  Temperature  vs  Radius  for  Various  Times  13 

4.  Variation  of  Bore  Convection  Coefficient  14 

5.  Variation  of  Conductivity  with  Temperature  15 

6.  Variation  of  Conductivity  with  Temperature  16 

7.  Variation  of  Specific  Heat  with  Temperature  17 

8.  Radial  Stress  Distribution  18 

9.  Hoop  Stress  Distribution  19 


I.  INTRODUCTION.  There  are  several  techniques  available  for 
quenching  metals.  The  object  is  to  develop  some  desired  micro- 
structure in  the  material.  In  the  Watervliet  Arsenal's  rotary 
forge  facility,  forged  cylindrical  tubes  are  heated  to  an  austen- 
itizing temperature  of  approximately  1550°F  and  then  quenched  to  form 
the  desired  martensitic  structure  in  the  material.  Both  external  and 
internal  diameter  quenches  are  utilized.  The  outside  diameter  is 
spray  quenched  with  four  water  jets  in  a diametral  plane  spraying 
water  onto  the  tube  while  the  tube  is  rotating.  There  are  several  of 
these  planes  located  along  the  axis.  The  bore  or  inside  diameter  is 
quenched  by  flushing  through  a nozzle  located  at  one  end  of  the  tube. 

While  the  facility  was  still  in  the  development  stage,  several 
tubes  (a  higher  incidence  than  normal)  developed  cracks  and  some  of 
these  were  interpreted  as  quench  cracks.  While  the  problem  was 
judged  to  be  metallurgical  in  nature  and  has  been  settled,  an  inter- 
est was  indicated  in  understanding  the  transient  temperatures  and 
stresses  involved  in  the  quenching  problem  and  this  led  to  the 
present  study. 

First,  the  transient  temperature  distribution  of  an  axially 
symmetric  hollow  cylinder  is  found.  Differences  along  the  axis  are 
assumed  to  be  minor  and  ignored.  The  thermal  properties  are  treated 
as  functions  of  temperature  rendering  the  equation  for  heat  conduction 
as  nonlinear.  The  finite  difference  method  is  used  to  solve  the 
temperature  problem.  The  Crank-Nicolson  equation  which  is  implicitly 
stable  is  used. 

In  the  present  study,  the  stresses  due  to  the  temperature  distri- 
bution and  the  martensite  transformation  were  computed,  assuming  the 
problem  was  elastic  and  linear.  From  the  computed  stresses,  it  was 
obvious  that  some  plastic  deformation  must  occur  and  that  an  elastic- 
plastic  analysis  was  required.  This  work  will  be  performed  in  a 
future  study. 

The  stresses  due  to  the  transformation  are  assumed  to  be  strictly 
due  to  a change  in  volume.  As  the  steel  transforms  from  the  austen- 
itic structure  to  a martensitic  structure,  a volume  increase  of  3%- 
4%  occurs  in  the  transformed  material.  This  volume  increase  gives 
rise  to  transformation  stresses. 
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The  thermal  problem  and  stress  problem  are  treated  as  being 
uncoupled.  The  heat  generated  from  the  transformation  is  small  and 
will  have  negligible  effect  on  the  temperature  distribution  during 
the  described  quenching  procedure.  The  transformation  begins  when 
the  temperature  in  the  material  reaches  Ms,  the  martensite  start 
temperature,  and  is  completed  when  the  material  is  past  the  Mf,  or 
martensite  finish  temperature.  Another  technique  used  in  quenching 
is  to  quench  the  material  to  Ms  and  slowly  allow  the  transformation 
to  take  place.  In  this  case,  the  heat  generated  during  the  trans- 
formation, might  be  significant  and  the  coupled  problem  might  need 
to  be  considered. 

II.  PROBLEM  STATEMENT.  The  partial  differential  equation  for 
temperatures  in  a hollow  cylinder  is 


1 3 rv  r ^u-»  r \ r \ (1) 

(k(u)r  " c(u)p(u)  at 

where  r represents  the  distance  along  a radius,  u the  temperature, 
and  t the  time.  The  thermal  conductivity,  specific  heat  and  density 
are  represented  by  k,  c and  p respectively.  These  properties  are 
assumed  to  be  functions  of  the  temperature.  Axial  symmetry  is 
assumed  and  any  effects  along  the  axis  are  ignored. 

The  initial  condition  is  given  by 

u(r,o)  = UD  (2) 

where  U0  would  represent  the  high  austenitizing  temperature.  The 
boundary  conditions  for  the  problem  described  would  be  of  the  con- 
vection type.  However,  to  allow  some  flexibility  in  the  program, 
they  were  written  in  the  following  form 

3u  , . 

37  - hlu  " -*i  at  r = a 

(3) 

3u  , 

37  " h2u  = -g2  at  r = b 

where  r = a specifies  the  inside  radius  and  r = b the  outside  radius 
of  the  cylinder.  The  values  h^  and  g^  can  be  varied  at  either  sur- 
face so  that  various  boundary  conditions  can  be  specified.  For 
example,  if  gj  = 0 at  r = a,  then  h^  is  the  Nusselt  number  and  the 
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convection  boundary  condition  is  indicated  at  r = a.  If  and  g2 
are  very  large  but  the  ratio  g2^2  = ^b  'then  the  temperature  Uj,  is 
specified  at  r = b. 


Since  the  thermal  properties  of  the  material  must  be  considered 
as  functions  of  temperature,  the  partial  differential  equation  (1) 
is  nonlinear  and  numerical  techniques  are  needed  to  solve  the  problem. 
An  implicit  scheme  based  on  the  Crank-Nicolson  equation  was  used  in 
writing  the  finite  difference  scheme  for  the  temperatures. 

III.  FINITE  DIFFERENCE  EQUATIONS.  The  Crank-Nicolson  repres- 
entation of  Eq.  (1)  is 


a+iAr 


(a+(i-%)Ar)Ar 


k • . , 6 (u.  , , + u.  , ) 

l+^.n+Ji  r i+h,n+l  i+h,n 


a+  (i-l)Ar  , P , 

— 2“  ki-%,n+Js6r(ui-%,n+l 


U, 


“i-ls.n’1  * 


c i ,n+'hf *i  ,n+h 


i , n+ 1 


- u • 


l.n 


At 


(4) 


where  i is  the  ith  node 
n is  the  time  step 
At  is  the  time  increment 
Ar  is  the  space  increment 


and 


6 u.  i 
r i+%, n 


ui+l,n  ~ ui,n 

Ar 


k = k 1 Ui.n+^ 

Ki+h,n+h  2 


(5) 


The  finite  difference  equation  (4)  is  written  about  the  point  r^, 
tn+%*  temperature  and  its  spacial  derivatives  can  be  written 

without  requiring  their  values  at  n+h,  then  the  equations  become 
linear  [1].  This  is  accomplished  by  arithmetic  averaging  the  finite 
difference  analogues  at  the  points  r^,  tn  and  r^,  tn+]_,  and  the 
resulting  analogue  is  the  average  of  the  forward  and  backward  ana- 
logues. To  solve  these  equations  for  the  temperatures,  it  is  required 


1.  von  Rosenberg,  D.  U.,  Methods  for  the  Numerical  Solution  of 

Differential  Equations,  American  Elsevier  Publishing  Co.,  Inc., 
New  York,  1969. 
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to  know  the  properties  at  the  n +h  time  step.  This  will  be  shown 
later. 

In  writing  the  boundary  conditions  as  Equations  (3) , it  is 
necessary  to  locate  the  nodes  as  shown  in  Figure  1.  There  are  no 
nodes  located  on  the  boundary  about  which  finite  difference  equations 
are  written.  The  boundary  conditions  are  used  to  eliminate  the  Oth 
and  R+lst  nodes  from  the  equations.  The  temperatures  on  the  boundary 
are  found  from  extrapolation  or  through  the  use  of  the  boundary 
conditions  after  the  spatial  temperature  distribution  is  found  at 
that  time  step. 

The  values  of  the  thermophysical  properties  at  the  h time  step 
can  be  found  through  various  projection  methods  [1] . The  one  chosen 
is  the  centered  Taylor  series  projection.  A set  of  equations  similar 
to  Equation  (4)  are  written  between  the  n and  n+h  time  step.  Thus 
the  values  of  the  properties  would  thus  be  required  at  the  n+h  time 
level.  The  technique  allows  the  computation  to  take  place  using 
properties  evaluated  at  the  known  nth  time  level.  Under  those 
conditions,  the  equations  are  linear,  the  properties  known  and  the 
temperatures  at  the  n+h  time  step  can  be  found.  Knowing  this,  new 
property  values  can  be  found  and  the  equations  solved  for  the  temper- 
atures at  the  n+1  time  step. 

An  alternate  technique  which  still  arrives  at  the  equivalent 
of  Equation  (4)  was  used  for  the  centered  Taylor  series  projection. 
Equation  (1)  is  rewritten  in  following  form. 


k(u) 


O’ 


kfu)  - ..  , , 8u 

T-  = P(U)C(U)  — 
r 8r  dt 


(6) 


Using  the  Crank-Nicolson  finite  difference  analogue  about  the  ith 
node  and  n+h  time  step,  Equation  (6)  becomes,  after  some  rearranging, 


^ r,  . . ...  . rui+l ,_n+h  ui-l,n+^.,  Ar 

2 (Ar)2  k Ui,n+Jp  +k  (ui,n+%^  2Ar  J 2^ui+l,n+l  ui+l,n^ 

1 n , s iif  •>  rui+l>n+%  - ui~l  ,n+%1  Ar-.  , . „ 

2 f Ar~> 2 l-k^ui,n+^  k ^ui,n+%^  2Ar  2 i_1»n+1  i_1»n^ 

(7) 


2(Ar) 


1.  von  Rosenberg,  D.  U. , Methods  for  the  Numerical  Solution  of 

Differential  Equations,  American  Elsevier  Publishing  Co.,  Inc., 
New  York,  1969. 


4 


(7) 


(Ar) 


^Tk(ui,n+^(ui,n+l  + ui,n^  + 


+ ui+l,n+l  * ui+l,n  ~ ui-l,n+l  ~ ui-l,n, 

2 (a+(i-%)Ar) 


2Ar 


= cCui,n+^PCui,n+^ 


ui,n+l  ~ ui,n 

At 


However,  the  coefficients  of  the  first  2 terms  can  be  viewed  as 
truncated  Taylor  series  for  kfui+j^n+j^) 

k fu  i s kfu  i + At  i. • / •>  ,-ui+l>n+%  ui~l>n+^-> 

^Ui+}5,n+V  ^ui,n+*r  2 ^Ui,n+V*-  2Ar  ^ 


C8) 


Rewriting  Eq.  (7) 

[kK+Js,n+^Hui+1>n+1  + ui+1>n)  + 


+ 2 (Ar) 2 tk(ui-Jj,n+Jp^ui-l,n+l  + ui-l,n^ 


2 (Ar)z  ^k^ui+J2>n+}P  + kfui-%,n+^^ui,n+l  + ui,n^  + 


C9) 


+ k(ui+^,n+^)  ^uj+i,n+l  * ui+l,n^  _ k(ui-%,n+^)  ^.ui-l,n+l  * ui-l,nj 
4(a+(i-%)Ar)  2Ar  4(a+(i-^)Ar)  2Ar 


^■c^ui+h,n+h^  + c^ui-%,n+Jp^P^ui+%,n+Jp  + P(ui-J5,n+J$)] 


[Ui>n+1  ~ ui,m  _ Q 
At 
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where 


kK,W  = %[k(ui+j£  + kCUi.j^n+jj)] 


(10) 


Equation  (9)  is  now  rewritten  for  evaluation  of  temperatures  at  the 
n+h  time  interval  using  for  the  coefficients  their  known  values  at 
time  step  n 


_1 

2 (Ar )1  2 


[1  + 


Ar 


4 (a+ (i-^)Ar) 


] 


k(ui+h,n)(ui+l,n+h  + ui+l ,n^ 


1 

+ 2 (Ar)2 


[1  - 


Ar 


4(a+(i-%)Ar) 


■]  k(ui-%,n)  (ui-l,n+%  + ui-l,n)  2 (Ar) z 


[k(ui+%,n)  + k(ui-^n^Hui>n+%  + ui>n)  - %[c(ui+%>n)  + c(Ui^>n)] 
[P(ui+Js,n)  + P(ui-%,n)]  — }%1-  — = 0 (11) 


At 


This  yields  the  temperature  distribution  at  n+h  using  coefficients 
evaluated  from  the  temperature  distribution  at  time  n.  The  temper- 
atures from  Eq.  (11)  are  then  used  to  evaluate  the  coefficients  for 
use  in  Eq.  (4). 

The  finite  difference  equations  were  tridiagonal  and  were  solved 
using  the  Thomas  algorithm  [1] . Briefly,  the  form  of  this  algorithm 
is 


aiui-i  + Vi  + Ciui+1  = di 

al  = 0,  cR  = 0 


l<i<R 


(12) 


where  the  terms  on  the  left  hand  side  are  at  the  n+1  time  step  and  on 
the  right  hand  side  at  n time  step.  For  all  i,  the  quantities 


3i  = bi  - 


aici-l 

h-i 


(13) 


1.  von  Rosenberg,  D.  U.  Methods  for  the  Numerical  Solution  of 

Differential  Equations,  American  Elsevier  Publishing  Co.,  Inc., 
New  York,  1969. 
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as  computed  and  then  back  substitution  is  used  to  find  the  temper- 
atures from 


UR  ~ Yr 
ui  = Yi 


ciui+l 

h 


(14) 


Computation  times  are  rapid. 

IV.  THERMAL  AND  TRANSFORMATION  STRESSES.  The  quenching  process 
gives  rise  to  thermal  stresses  due  to  the  large  thermal  gradients 
that  exist.  Areas  near  the  boundary  are  cooler  than  interior  points. 
The  boundary  would  like  to  contract  but  is  partially  prevented  from 
doing  so  because  of  the  interior,  hence  tensile  stresses  are  set  up 
near  the  boundaries  while  the  interior  is  in  compression.  The 
thermal  stresses  in  an  axially  symmetric  hollow  cylinder  subject  to 
a non-uniform  temperature  distribution  are  given  by  [2] . 

Fn  _2  2 P ,r 

p_  = rr  [~o — t / Pu(P)dP  " / pu(p) dp] 

r b2-a2  a a 

(15) 

p 2 2 b r 

Pfl  = [— — - ■ / Pu(p)dp  - / pu(p)dp  - r2u(r)] 

r b2-a2  a a 


where  Pr the  radial  stresses 

Pq the  tangential  stresses 

E Young's  Modulus 

a thermal  expansion  coefficient 


u(p) .. .radial  temperature  distribution 

The  stresses  due  to  the  transformation  are  found  using  similar 
equations  since  these  stresses  are  due  mainly  to  a volume  increase 
in  the  transformed  material.  The  difference  between  the  two  calcu- 
lations is  that  the  transformation  does  not  occur  across  the  thick- 
ness simultaneously  but  progresses  across  based  on  the  temperature 
in  the  cross  sections.  No  transformation  stresses  exist  when  the 
temperatures  are  all  above  that  temperature  (M_)  when  the  trans- 
formation begins  or  below  that  temperature  (Mf)  for  which  the 
transformation  ends.  Between  these  two  temperatures,  a linear 
change  in  volume  is  assumed.  The  change  in  volume,  about  4%  if 
the  transformation  is  complete,  is  assumed  to  be  isotropic  so  that 


2.  Boley,  B.  A.  and  Weiner,  J.  H. , Theory  of  Thermal  Stresses, 
John  Wiley  and  Sons,  Inc.,  1960. 
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it  translates  to  one-third  of  the  volume  change  for  a linear  change 
during  the  transformation.  Stresses  are  computed  in  a manner 
similar  to  thermal  stresses. 


(16) 


where 


I. 


j-  = 0 if  u(r)  >_  Ms 

A / 


(17) 


The  Ms  temperature  was  taken  to  be  350°F  and  the  Mf  temperature  150°F 
for  the  computations.  Figure  2 shows  some  temperature  distributions 
which  can  arise.  In  the  upper  figure,  no  transformation  has  taken 
place,  hence  the  transformation  stresses  are  zero.  In  the  lower 
figure,  the  transformation  is  occurring  from  both  the  inside  and 
outside  radius.  In  sections  indicated  by  I,  corresponding  to  Equa- 
tions (17),  the  transformation  has  not  begun,  in  sections  II,  the 
transformation  is  progressing  and  in  sections  III,  the  transformation 
is  complete. 

As  stated  above,  for  the  present  study  the  thermal  and  trans- 
formation stresses  were  assumed  to  be  elastic.  In  the  following, 
the  results  indicate  that  stresses  are  too  large  for  this  assumption 
to  be  valid.  References  [3]  and  [4]  treat  similar  problems  using 
elastic-plastic  analysis. 


3.  Weiner,  J.  H.,  and  Huddleston,  J.  V.,  Transient  and  Residual 
Stresses  in  Heat-Treated  Cylinders,  Journal  of  Applied  Mechanics, 
March  1959. 

4.  Landau,  H.  G. , and  Zwicky,  E.  E.  Jr.,  Transient  and  Residual 
Thermal  Stresses  in  an  Elastic-Plastic  Cylinder,  Journal  of 
Applied  Mechanics,  September  1960. 
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V.  RESULTS  AND  DISCUSSION.  Figures  3 through  7 show  some 
resulting  temperature  distributions  under  various  conditions . 

Figure  3 shows  the  temperature  distribution  across  the  wall 
thickness  for  various  times.  The  boundary  conditions  are  convective 
and  h2  represents  a value  near  that  suggested  by  the  manufacturer 
of  the  quenching  facility  for  the  coefficient  on  the  outside  diameter. 
On  the  inside  diameter,  the  value  of  h^  was  said  to  be  lower  than 
that  of  h2«  The  radius  is  in  inches  and  the  temperatures  are  in  °F. 

The  ambient  temperature  was  assumed  zero.  Figure  4 shows  the  effect 
of  variations  in  the  convection  coefficient  on  the  inside  diameter. 

The  results  are  shown  for  only  one  time  step.  Since  the  temperature 
at  that  time  does  not  change  much  under  the  different  boundary  condi- 
tions, small  differences  in  the  convection  parameter  on  the  inside 
diameter  will  have  little  effect  on  the  transformation. 

Figure  5 shows  the  effect  when  the  thermal  conductivity  is 
allowed  to  vary  with  temperature.  For  the  same  time  step,  three 
curves  are  shown.  The  conductivity  is  allowed  to  be  an  increasing, 
decreasing  or  constant  function  of  temperature.  Finding  real  data 
to  use  in  the  program  is  difficult.  Figure  6 shows  the  temperature 
distribution  for  a bilinear  thermal  conductivity  curve  based  on  one 
for  4130  steel.  These  properties  are  usually  determined  experimentally 
under  equilibrium  conditions.  Since  the  structure  of  the  material 
is  changing  under  rapid  cooling  and  since  equilibrium  does  not  exist, 
the  properties  which  should  be  used  are  those  determined  under  the 
same  conditions  as  the  quench.  This  can  be  described  best  by  looking 
at  the  specific  heat.  Again  for  4130  steel,  a spike  increase  in  the 
value  of  the  specific  heat  occurs  between  1200°F  and  1500°F.  This 
occurs  during  heating  and  is  due  to  the  austenitizing  of  the  material. 
Figure  7 shows  the  temperature  distribution  throughout  the  tube  wall, 
allowing  the  specific  heat  to  be  a function  of  the  temperature  but 
ignoring  the  spike  mentioned  above.  Under  the  quench  conditions,  the 
spike  would  occur  during  the  martensite  transformation  and  be  of 
different  shape. 

Figures  8 and  9 show  a sample  of  the  radial  and  circumferential 
or  hoop  stress,  respectively,  for  a specific  time.  They  are  taken 
from  one  of  the  previously  cited  cases.  The  insert  shows  qualita- 
tively the  temperature  distribution  at  the  time  indicating  that  the 
transformation  has  begun  at  both  the  inside  and  outside  radius.  Each 
figure  shows  the  thermal,  transformation  and  sum  of  the  stresses. 

Large  compressive  hoop  stresses  indicate  the  strong  possibility  of 
plastic  deformation.  Since  the  transformation  occurs  at  the  lower 

temperatures,  the  thermal  gradients  are  smaller  and  the  thermal 
stresses  lower  than  their  values  earlier  in  the  quenching  cycle. 

The  thermal  stresses,  however,  can  be  a significant  part  of  the 
total  stresses,  especially  early  in  the  transformation,  and  should 
not  be  neglected  in  an  elastic-plastic  analysis. 
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FIGURE  2.  TYPICAL  TEMPERATURE  DISTRIBUTIONS 
NEAR  TRANSFORMATION  TEMPERATURE 
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FIGURE  4.  VARIATION  OF  BORE  CONVECTION  COEFFICIENT 


14 


RADIUS,  inch 
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